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We theoretically study electronic transport through a layer of quantum dots connecting two metal- 
lic leads. By the inclusion of an inductor in series with the junction, we show that steady electronic 
transport in such a system may be unstable with respect to temporal oscillations caused by an 
interplay between the Coulomb blockade of tunneling and spin accumulation in the dots. When this 
instability occurs, a new stable regime is reached, where the average spin and charge in the dots 
oscillate periodically in time. The frequency of these oscillations is typically of the order of 1GHz 
for realistic values of the junction parameters. 



I. INTRODUCTION 



the junction. 



During the last decade, spin-polarized electronic trans- 
port through quantum dots (QD) connecting two metal- 
lic ferromagnetic leads has been the subject of very in- 
tensive both theoretical and experimental research with 
numerous applications in spin based devices such as spin 
valves, spin filters, spin diodes, etc. 1- ''' Ferromagnetic or- 
dering in the leads, which causes the rates for tunneling 
on-to and out-from the dot to depend on the electron 
spin, results in the accumulation of spin in the dots 4 ' 5 
in addition to charge. This spin accumulation affects 
the transport properties of the system*''' and provides 
a way to control the spin polarization of the current by 
bias or gate voltages. Moreover, it was found that the 
Coulomb blockade (CB) phenomenon significantly affects 
the DC-current 8-11 and the shot noise 12 in such systems, 
which has provided new opportunities for their electrical 
manipulation 13 . 

In most previous papers, time- independent effects 
caused by the spin and charge accumulations on the dots 
between ferromagnetic leads have been investigated. In 
this paper we focus on self induced time-dependent phe- 
nomena which may arise in such systems under a static 
bias with an inductor included in series with the junc- 
tion. We will show that under certain conditions, time 
independent electronic transport across a layer of quan- 
tum dots placed between normal and magnetic leads be- 
comes unstable, which results in oscillations in time of 
spin and charge accumulations in the layer of dots. One 
may assign this instability to the existence of a negative 
differential conductance (NDC), arising due to an inter- 
play between the Coulomb blockade and the spin block- 
ade phenomena. The occurrence of an NDC in systems 
similar to the one considered by us was reported in 14-18 . 
However, we will show that in our case, an instability of 
the time-independent regime of the charge and spin flow 
may arise even at a positive differential conductance of 



II. MODEL 

We consider a layer of identical quantum dots, with 
two spin- dependent states, connecting normal (N) and 
ferromagnetic (F) metallic leads. All dots in the layer are 
supposed to be placed at the same distance with respect 
to the leads (see Fig. 1). The system under consideration 
is described by the Hamiltonian 



H = J2 H l t ead + H QD + J2 H i 



tunnel 



(i) 



i—N,F 

where the partial Hamiltonians 



i=N.F 



PI 



H QD 

tunnel 



EEr at ar -, 
k,a,i fc)<r»* 5 

k.rr 

e cr<4-,n C c,n + ^ c },n C t,^ C l,n C ii^' (^) 

a. n 

T ^y^[ a r Ca.n+h.C.] 



describe electrons in the leads, in the QDs, and tunneling 
coupling between QDs and leads, respectively. Here at 

creates an electron with wave vector k and spin a = J [ 1 J, 
in the corresponding lead i = TV, F (i is the lead index) ; 
Er N = e(k) and E% F = e(k) — I a where e(k) is elec- 
tron kinetic energy and 1^ = — = I is the ferromag- 
netic exchange energy; c\ n creates an electron with spin 
a and energy e CT in the n-th dot (n = 1,2, J\f) counted 
from the identical chemical potentials of electrons in the 
metallic leads fj,o] U is the Coulomb interaction energy 
due to the double occupancy of the QD level by electrons 
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FIG. 1: Charge and spin transport through a layer of quantum 
dots sandwiched between a normal (N) and a ferromagnetic 
(F) lead subjected to an external magnetic field B is studied. 
The junction is symmetrically biased by the voltage 2Vb- 



with opposite spins. We consider the case when the inter- 
dot distance is much larger than the distance between 
the metallic leads allowing us to neglect the inter-dot 
Coulomb interactions. We set /irj as an origin for mea- 
suring all energies. The difference between spin "up" and 
"down" energy levels Ae = e^-— ej, in QD can be controlled, 
e.g. by the Zeeman splitting induced by the applied ex- 
ternal magnetic field B which we will fix, for the sake of 
definiteness, in "down" direction (Ae « 10~ 2 -r- 10 _1 meV 
in the magnetic field 0.1 -r- IT 19 ). The external magnetic 
field and the F-lead magnetization are taken parallel in 
this model, so no spin precession effects, as investigated, 
e.g., in Ref. 20, are present. We will consider both the 
possible orientations of F-lead magnetization: along the 
magnetic field direction (P) and opposite to it (A). Fur- 
thermore, due to the exponential sensitivity of the tun- 
neling matrix elements to the geometrical position of 
the QD with respect to the leads, physically interesting 
limiting cases can be achieved. We have studied two such 
cases: the N-junction, in which the QDs are closer to the 
N-lead (tn 3> Tf), and the F-j unction where they are 
closer to the F-lead (tf 3> tn). 

To study charge and spin transfer between leads we 
will consider, for the sake of simplicity, a symmetric volt- 
age biasing of the junction (see Fig. 1). In this case the 
bias voltage does not affect the position of the QD lev- 
els e a , but it shifts the chemical potentials in the leads 
(see Fig. 2). In this paper we will restrict our study to 
the case when e^ > 0. The magnitude of the splitting 
Ae = e^ — ej, is assumed to be substantially larger than 
the intrinsic level width and ksT <C Ae (T is tempera- 
ture) thus providing well controlled separation of the spin 
states. Furthermore we will consider Coulomb blockade 
regime when temperature and bias voltage are much less 
then the double charging energy U. The conditions above 
determine the lack of charge and spin transfer if the ab- 
solute value of bias voltage |V&| is less then V a = |e CT /e| 
(e is the electron charge). When \Vb\ exceeds Vj., the J.- 
spin polarized tunneling arise manifesting itself as a step 
in the IVC (current- voltage characteristic) at |V;,| = V\,. 
By further increase of the bias voltage an additional step 
appears at \Vb\ = Vf when f-spin polarized tunneling is 
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FIG. 2: Electronic energy scheme (P-configuration) showing 
the nonmagnetic band (N) and the bands for \ (minority) and 
4- (majority) spins in the ferromagnet (F). The latter are split 
by 21, where I is the ferromagnetic exchange energy. The 
applied bias voltage 2V& shifts the electron energy by ±eVb in 
N/F-lead respectively from the chemical potential fio. Here 
e± and e-f are the spin-dependent energy levels in the dots. 



triggered. However, due to the Coulomb blockade, the 
presence of an electron in one state effectively blocks the 
current through the other one, thus affecting the IVC. We 
will show that a negative differential conductance can be 
achieved in the vicinity of Vf if the tunneling rate of the 
electron transfer between the f-state in the dot and the 
drain electrode is low enough compared to the one for 
the ^-state (spin blockade). 17 

To analyze quantitatively kinetic properties of the sys- 
tem under consideration we will use rate equations for 
the probabilities P|-(n to find electron with spin j" (I) 
on the dot. They can be derived from the generalized 
master equation for the density matrix by adopting the 
Markovian approximation in the limit of weak tunneling 
(see Appendix). 1 ' 1,21 Double occupation of the dots is pro- 
hibited by the Coulomb blockade thus leaving only two 
independent components in the rate equations (the prob- 
ability to find the dot unoccupied is Pq = 1 — Pj- — P|_). 
The combinations P c = Pj- + Pj, and P s = P|- — Pj, then 
describe charge and spin accumulation in the layer, re- 
spectively. The rate equations describing time evolution 
of the average dot populations can be presented in the 
form 

1 dP, 



Tm dt 



° r = -a + la )p a + (r N (v) + imv)) (i - p_ CT ). 

(3) 

Here Tjv = 2ir\TN\ 2 gN 1 gjy is the density of electronic 
states in the normal lead, which we assume to be energy 
independent; 
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where g a F is the spin dependent density of electronic 
states in the ferromagnetic lead; 
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(a) 



N-junction 



MV) = f F *(-V) = f(e a -eV), 

where 2V is the total voltage drop across the junction. 

The average current per one dot through the N— >QD 
junction is given by the expression 



j(V) = eT N J2 [-Pa + JnAV) (1 - P-a)] ■ (6) 



The geometrical asymmetry parameter £ = \tf/tn\ 2 
classifies 2 types of junctions: N-junction for £ <C 1 and 
F-junction for £ 3> 1. 



III. THE DC TRANSPORT 

Solving the system (3) at V = Vb = const we obtain 
the time independent occupation probabilities 



(1 + 7-cr)(/jV<7 + fFala) ~ ]\U N a' + Ifo'Io') 
(V 

n(i+7 CT ')-n(jw+iW'7<r') 



(7) 

Inserting expression (7) into (6), one gets the DC- 
current per dot, jo(Vb), through the junction. The IVC 
together with corresponding charge and spin accumula- 
tion probabilities are presented in Fig. 3. 

From the expression (6) and Fig. 3 it is evident that 
the IVC is not symmetric with respect to a change of 
sign of the bias voltage. In our further investigation we 
will focus our attention on the phenomena occurring 
by activation of the spin-up channel of the electron 
transfer. Therefore, we will consider the bias voltage 
values \Vb\ ps V\. Taking into account the fact that 
\eVb\ — ei « Ae 3> fcsT, one can distinguish two cases: 

a) NF-biasing (eV b > 0; f Fa {V b ) « 0, f Nl {V b ) « 1) 

b) FN-biasing (eV b < 0; f N<r (V b ) « 0, f n (V b ) « 1), 

in which we obtain an approximate expression for 
DC-current, accurate up to the exponentially small 
correction exp(— Ae/ksT) <C 1. 

a) NF-biasing 



jo(Vb) ~ er7v7l 



1 +7t~/A f t( 1 ~7t) 
(1 +7f)(l +74-) - fm 



(8) 



From Eq. (8) it is clearly seen that the voltage de- 
pendence of the stationary IVC is entirely driven by the 
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FIG. 3: Single-dot IVC with corresponding on-dot charge (P c ) 
and spin (P s ) accumulation probabilities for: (a) an N-type 
junction with j-f = 0.085, 7^ = 0.115 (£ = 0.1), and (b) an 
F-type junction with 7 t = 8.5, 74. = 11.5 (£ = 10). The 
current is normalized as jo — jo/eFN- We take e-f = 2ej_, 
temperature ksT/Ae — 0.05 and introduce the dimensionless 
voltages V — eV/Ae; V a = e CT /Ae. One can see that the IVC 
is not symmetric with respect to voltage biasing Vb — > —Vb- 
Also, in NF-biasing, the N-type junction produces a NDC 
around Vf — 2 in contrast to the F-type junction. For FN- 
biasing the NDC is absent in both cases. 



Fermi function ,/Wt(14) which exhibits a step- up like be- 
haviour at |Vb| = V^. Differentiating Eq. (8) with respect 
to voltage we find the differential conductance 



djo 
dV b 

where 



2 r ft 7j(l +7t) q, S 

[(1 +7 t )(l + 7;) - f m \ 



^(7T>7i) 
/jvt 



7T - 74 
df 



de 



e t -eV b 



7t74. 
< 0. 



(9) 



(10) 



The corresponding current "jump" Ajo = jo(— V\ 
SV) - jo(-V t - SV), where Ae > \eSV\ > k B T, is 



Aj = eT 



N - 



74. 



-tf(7t,74-)> (11) 



(1 +7i)(7t +74 +7t7i) 
in which $(7-)-, 74.) determines the sign of the "jump". 

From Eq. (9) it can be seen that the IVC ex- 
hibits NDC, djo/dVb < 0, if 7-)- and 74 satisfy condition 
i?(7t,7j.) < shown graphically in Fig. 4. To lowest 
order in the small parameter formed by the ratio of the 
polarization parameter and the Fermi energy, I/ef <C 1, 
this inequality may be written as 
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FIG. 4: (a) An IVC-section with a NDC around V b = Vf ap- 
pears for parameters 74- and 74. within the range correspond- 
ing to the shaded area based on the condition i?(7t, 74.) < 0. 
(b) The same condition expressed in terms of the normalized 
densities of states g a = gp/gN for various values of the ge- 
ometrical asymmetry parameter £ characterizing N-j unctions 
(£ < 1) and F-junctions (£ > 1). The diagonal line divides 
the region of parameters for the P-configuration (g^ > g^) 
and the A-configuration (<L < <?^). Evidently, NDC exists in 
P-configuration only. 



Therefore, the smaller is £, the smaller ferromagnetic 
polarization of F-lead is required for NDC. Obviously, 
NDC is present only in P-configuration. As it is seen 
from Eqs. (5) and (14) the differential conductance (9) is 



proportional to T 
decreases at \V b J 
"jump" in IVC is - 

b) FN-biasing 



1 at V b = 
Vf\ > 
k B T. 



— V4. Since it exponentially 
k B T/\e\, the width of the 
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(14) 



where 



/f4 — 



< 0. 



(15) 



One can see that NDC does not appear in the FN-biasing 
of the junction. 

The asymmetry of IVC with respect to the direction of 
bias voltage for junctions with magnetic leads is known 
and utilized e.g. for spin filtering, spin diodes etc. 2 ' 3 ' 18 
In the system under consideration NDC also occurs only 
at the regime of an NF-biasing. Physics of this phe- 
nomenon, observed experimentally 1 ', is most transpar- 
ent in the case of a strongly asymmetric N-P-junction 




FIG. 5: Circuit scheme used for the instability investigation: 
C is the intrinsic junction capacitance, C is the inductance of 
the circuit, J is the current flowing through the inductor. 



(7t < C7|< 1) at low temperatures. Indeed, under such 
conditions and at ej, < eVb < e-j- spin up states in the 
dots are not populated (P-f = 0), while those for spin 
down are almost completely populated (Pj, ~ 1), giving 
the average current jo ~ eTj^y^. Then at eVj, > e-f the 
up spin states start to contribute and the electrons are 
mostly trapped in these states (P-f ~ 1) due to a very long 
escape time ~ ^Fat) -1 that blocks the current through 
the spin down states (Pj. — Tf/T4- ^ 1) by the Coulomb 
blockade effect. Consequently, for eVb > e-f the average 
current jo ~ erjv(7fP-|- + 74. ^4) — 2eFAr7^ is less than 
the average current at eVb < £4- 



IV. THE DYNAMICAL INSTABILITY OF 
CHARGE-SPIN ACCUMULATIONS 

In this section we will only consider an iV-junction 
(la < 1) in the regime of NF-biasing which is the most 
interesting one. Really, a voltage biased circuit with a 
nonlinear resistor, R(V), providing the section with NDC 
in the IVC, is usually considered as a prerequisite for dy- 
namical instability of the system. To investigate stability 
of charge/spin accumulation in our system we introduce 
an inductor with inductance C in series with the junction 
(see Fig. 5). We also take into account the capacitance 
of the junction C w NC\qd, where M is the number 
of quantum dots inside the junction, C\qd = sl 2 /d is 
the average capacitance per one dot (e is the dielectric 
constant of the layer material, I is the average distance 
between dots, d is the distance between the leads). 

The tunneling and electrical processes in this system 
are governed by the system of four differential equations: 



1 dP a 
Tn dt 



(1 + la )P a - f N a(V)(l - P- a ) = 

dj 



C—+2V = 2V b (16) 



2C^ = J. 
dt 



A/j(V,P CT 

Here J is the total current flowing through the indue- 
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tor, and the current per one dot j(V, P a ) is given by Eq. 
(6). The set of differential equations (16) always has the 
time-independent solution 

P a = P°\ V = V b , J = Mj(V b ,P?) (17) 

with time-independent probabilities P§* given by expres- 
sion (7). We analyze the stability of this solution by 
linearizing the set of equations in terms of small devia- 
tions ~ exp(Ai) from the time- independent solution (17), 
considering the bias voltage V b ~ Vf and the tempera- 
ture k B T < Ae being low enough to fulfill fp(V b ) « 0, 

fh{V b ) « 1. Thus transfo rming system (16) into a set of 
algebraic equations, we obtain the characteristic equation 
of the fourth order in the Lyapunov exponent A, i.e. 



+ a 3 X 3 + a^A 2 + OiA + ao = 0, 



(18) 



with the coefficients 

Oi = T N A 2 -r i~JV ei 

a 2 = l + £AfeT 2 N A 3 ( lf + ll ) 
a 3 = CNeT N A 3 + CCT N A 2 

04 = CC 



£Afer%A 3 $ 



CCT 2 N A l 



(19) 



where 



Ai 
A 2 



(1 

24 



h7t)(l + 
7t + 74 



74.) _ fm 



(20) 



74-( 1 + 7t) 



A't 



(l + 7 t )(l + 74.) - f Nf 



The analysis of the characteristic polynomial in the 
complex A-plane, namely counting the winding number 
of its phase as the variable encircles the Re(X) > half- 
plane, shows that two roots of the equation (18) are al- 
ways real and negative, while the other two have the real 
part changing its sign as C passes through a critical value 
C c . In order to find the critical values of the inductance 
at which the real part of two complex conjugated Lya- 
punov exponents is equal to 0, we insert A = iio (were ui 
is real) into the characteristic equation (18) and obtain 
a set of two equations 



d4(£)u; 4 — a 2 (C)uj 2 + ao 
a 3 (C)uj 2 - ai(£) 



(21) 



from which the critical inductance C c and the corre- 
sponding roots Ai 2 = ±iw c of the characteristic equation 
are found. 

When the tunneling rates, which control the aver- 
age electron populations of the dots, arc the fastest 



rates in the system, NDC leads to the well-known 
electro-dynamical instability of the time- independent 
current flow. It is interesting to note that in the 
system under consideration an instability of the time- 
independent regime arises even in the case when the 
RC-time t^c = \Rd\C, where the differential resistance 
Rd = (Afdjo/dVb)~ 1 is defined by Eq. (9), is the shortest 
time scale in the system: 



tRC <C r« , t RC <C 



(22) 



Taking into account realistic values of the junction 
parameters I ~ lOnm, d ~ lnm, e ~ 10~ n F/m and the 
operating temperature T ~ 50mK, one finds that the 
inequalities (22) are well fulfilled down to j a > 1CP 2 and 
inductances CAf^ > 10" 10 i7. 

Conditions (22) permit skipping of all terms containing 
C in system (16) and coefficients (19), i.e. putting there 
C = 0. In this case the set of equations (21) reduces to 
an algebraic quadratic equation for C with roots 



£(±) 



Rd 



^±^ 2 -4[(l + 7+ )2-l]tf 



T N [(1 + 7 t )(l + 7i ) - f m h^ 



(23) 



where $(7^, 74,) is defined by Eq. (10); 7+ = 7-f + 74, and 
^(7t> 74, V b ) = 1 - /jvt ~ 27 t - 

As it follows from Eq. (23), in the case of Rd < 
(and hence ■& < 0), one root is negative, and the other 
one is positive. Therefore, there is only one critical 
value of the inductance C c = at which the system 
looses its stability. On the other hand, easily seen from 
coefficients (19), in the absence of external inductance 
(C = 0) the characteristic equation (18) reduces to 
the second order with both solutions having real parts 
Re(X) < 0, indicating that the fixed point is stable for 
any choice of the other parameters. From there it follows 
that in the range of parameters in which the differential 
resistance Rd is negative (light grey area in Fig. 6), the 
time- independent solution (17) is stable if < C < C c , 
and, if the inductance exceeds the critical value L > C c , 
the system looses its stability. 

Following from Eq. (23), the presented system reveals 
one peculiar property: it may also loose its stability in 
the case of a positive differential resistance Rd > (that 
is d > 0). In the range of 7 CT parameters satisfying 



1) > 0; Lp 2 



> 



4 [(1 + 7t + 74) 2 - 1] ^ 



(24) 



(shown as a dark grey area in Fig. 6) both roots are 
real and positive and hence there are two critical values 
of the inductance C c = and £? c = One can 

find that the system is unstable for inductance laying in 
the interval £ c < £ < C° and stable otherwise. 
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FIG. 6: The instability condition in (7^ , 74_)-parameter space: 
the light grey area shows the range of parameters in which 
Rd < and the instability arises for C > C c . The dark 
grey area shows the range of parameters in which Rd > 
and the instability arises for C c < C < Here Vb — Vf, 
fc s T/Ae = 0.05. 



All the features mentioned above are most transparent 
in the limit of a highly asymmetric N-junction 7, C 1. 
In this case the critical values of the inductances C c and 
as well as the corresponding frequencies uj c = w(£ c ) 
and = w(£c)> are 



2 (l - fm? 

er 2 N A% T 7i(7T + 7i)^ 



UJ,, 



r N Vtf. (25) 



From expressions (25) it is clearly seen that in the limit 
7cr <C 1 the first critical value (£ c ; oj c ) does not depend 
on sign of differential resistance, while the second one 
(£*; w*) appear as soon as IVC attains positive differen- 
tial conductance ($ > 0). Therefore, the instability takes 
place both for Rd < and Rd > as soon as the induc- 
tance exceeds C c , but, as mentioned before, stability is 
established again as soon as C > if the differential 
resistance is positive. 

The above-mentioned instability is a Hopf bifurca- 
tion, resulting in the onset of spontaneous, non-linear, 
periodic in time, self-excited oscillations of current J(t), 
voltage drop V(t), average charge q(t) = eNP c {t) and 
average spin s(i) = (1/2)N P s (t) in the layer of dots. 
In the case (C — C c )/C c <C 1 the frequency of the 
oscillations w uj c . Analytical estimations and numerical 
calculations show that the critical inductance increases 
and the oscillations fade out with an increase of the 
temperature, disappearing at ksT ~ Ae. The numerical 
solutions of system (16) with C = 0, for an N-type junc- 
tion in P-configuration with the NDC, in circuit with 
inductance chosen slightly beyond the critical value, 
show charge/spin accumulation, as well as current, 




V-V n 



FIG. 7: (a) The limit cycle in (V, P c , P s )-space around the 
fixed point (0.5,0.905,-0.743). (b) The limit cycle in the 
(V, J)-space around the fixed point (0.5,0.102) at the IVC 
(gray curve). We use scales: V = eV/Ae, J = J/eTN-Sf, 
C — Ce 2 Y 2 N M I Ae, oJ = uj/Yn- Voltage is measured with 
respect to Vo = {Vf + V\)/2. Here the choice of parameters 
is: V b = F t , fcsT/Ae = 0.05, 7 t = 0.085, 7^= 0.115, £ = 0.1, 
C — 8.0 with critical values for instability £ c — 7.8712, ZZ? C = 
0.5457. This choice of parameters (P-configuration) leads to 
an IVC with NDC and the frequency of stable oscillations 
U= 0.5456. 



orbiting in time along the limit cycle as presented in 
Fig. 7. An example of limit cycle for A-configuration 
with positive differential conductance is shown in Fig. 
8. The results show critical inductances of the order of 
lmH/A^ and critical frequencies of the order of Tn, i-e. 
1GHz, independent of Af. Typical junction capacitance 
per dot C\qd ~ 10 _18 F introduces just small correction 
to critical inductance lowering it by ~ 0.2 percent, while 
critical frequency remains unchanged. 



V. CONCLUSION 

We considered the effect of Coulomb blockade correla- 
tions on the spin-dependent electronic transport across a 
layer of quantum dots connecting a normal and a mag- 
netic lead. It was shown that in such a system, under the 
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voltage biasing with an inductor added in series with the 
junction, an instability in the steady (time-independent) 
flows of charge and spin may arise. This instability de- 
velops into a new stable regime in which the average spin 
and charge accumulated in the dots oscillate periodically 
in time. The typical frequency of the oscillations is of 
the order of 1GHz for realistic junction parameters. In 
contrast to the standard electric instability of an RLC- 
circuit with a negative differential resistance, in the sys- 
tem under consideration spin accumulation in the dots 
results in an instability which occurs even in the case of 
an RL-circuit with a positive differential resistance. 
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VI. APPENDIX 

FIG. 8: The limit cycle appearing as the result of the insta- 
bility in the system with positive differential conductance at 
the stationary IVC (gray curve) in A-configuration, around a) 
the stationary points (0.5,0.981, —0.948) in (V,P C ,P S ) space, 
and b) around stationary point (0.5, 0.0206) in (V, J) space. 
Used scales and voltage origin are the same as in Fig. 7. The 
choice of parameters is: Vb = V^, kBT/Ae = 0.05, jf = 0.08, 
74. = 0.02, £ = 45.0. The critical values between which the 
instability takes place are £ c = 41.112; u; c = 0.5883 and 
Z° c = 257.301; uf c = 0.3272. The frequency of the stable 

oscillations uj = 0.5865. The rate equation 



Pi J l= N,F \ r l,L^QD Q _ T i.QD^L \ P l 



where eight T's are the rates of electron tun- 
neling between the lead ("L") and the dot 
("QD"), calculated using Fermi golden rule 
and written in compact form as Tl^iVi) = 
xN 2 / /" ( E iAP) + eVi - W) S (E it *(p) + eV t - e a ) dp, 
resulting in the expression 



rr l (V) = -\T l \ 2 p(e a + , l eV- f , ) 
xgi (et + e„ - ^(1 + ^)al + neV^j . (27) 

Here index i = N,F denotes process between QD and 
N/F lead, a = ±1 stands for spin f / J,, q = =pl for 
i = N/F, while 77 = L -> QD,QD -> L denotes the 
tunneling process from the lead to the QD and vice versa 
respectively. Tj is averaged, energy independent matrix 
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element of tunneling Hamiltonian (2) between QD and 
the lead i, gi(e) is electron density of states in lead i at 
given energy e, while f v (e) denotes Fermi function /(e) 
for L — > QD process, and 1 — /(e) for QD — > L pro- 
cess respectively. Since g(V) is slowly varying function 
along the voltage interval (Vj,,V^-), we approximate the 
tunneling rates as 



where Tf = (2ir / h)\Ti\ 2 gf is the voltage-independent 
"bare" tunneling rate, by approximating g? g(£i — ^(1 + 
<,i)al) for eV and e a — fi e^,/ (e<x — can be adjusted 
small using the gate voltage). 



r i ^(v)«/"( Cff + fte v- W) )rj 



(28) 
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